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We critically reexamine the bosonization-debosonization procedure for systems including certain 
types of localized features (although more general scenarios are possible). By focusing on the case of 
a tunneling junction out of equilibrium, we show that the conventional approach gives results that 
are not consistent with the exact solution of the problem even at the qualitative level. We identify 
inconsistencies that can adversely affect the results of all types of calculations. We subsequently 
show a way to avoid these and proceed consistently. The extended framework that we develop here 
should be widely applicable. 


I. INTRODUCTION 

It is well known that quantum physics gets richer and 
more peculiar as one considers reduced-dimensionality 
scenarios. Such scenarios are nowadays far from be¬ 
ing esoteric. From highly anisotropic and artificially 
layered materials, to nanostructures, to confined ultra¬ 
cold atomic gases, examples abound of what was once a 
playground for theorists but modern experimental tech¬ 
niques turn into a practical reality. The physics of one¬ 
dimensional systems provides an example in which strong 
quantum effects together with interactions and restricted 
kinematics modify the expectations we bring with us 
from our more familiar three-dimensional world [1, 2]. 
A case in point, the successful paradigm of the Landau 
Fermi liquid generically breaks down and gives rise to a 
new type of quantum fluid known as the Luttinger liquid 
[3]. A technical stepping stone on which the generality 
of this new paradigm rests is the technique known as 
bosonization. 

The term bosonization refers to the practical possibil¬ 
ity of describing the excitations of fermionic systems via 
a description based on bosonic degrees of freedom. The 
key observation is that for a fermionic one-dimensional 
system with strictly linear dispersion and no cutoff, the 
excitations at constant fermion number are particle-hole 
pairs that can be used to construct bosonic operators 
which completely capture the full excitation spectrum; 
such a view is known as the constructive approach [4]. 
The conceptual advantage of the constructive point of 
view is that it highlights the fact that bosonization is an 
exact correspondence between the two systems. There 
are also various complementary presentations based on 
the matching of correlators and known as the field- 
theoretic or the hydrodynamic approaches (these are, for 
example, more amenable to the conceptual description of 
the phenomenology of Luttinger liquids) [1, 2, 5]. 

Since the conditions of linear dispersion with large 
bandwidth and conserved particle number are all nat¬ 
ural approximations for systems at sufficiently low tem¬ 
peratures (much lower than the bandwidth and the en¬ 
ergy range of deviation from linearity), the applicabil¬ 
ity of bosonization is ubiquitous for all types of one¬ 


dimensional systems. Moreover, since within its appli¬ 
cability conditions the bosonization mapping is exact, it 
can equally well be used in both equilibrium and out-of¬ 
equilibrium situations. In particular, it provides a fertile 
ground for the study of transport phenomena in a variety 
of settings. 


A. Transport Problems and Bosonization 

Due to the versatility of the technique, a sizable frac¬ 
tion of the current theoretical studies of transport in 
one-dimensional settings rests on the use of bosoniza¬ 
tion. If the applied voltages are sufficiently low (in the 
same sense as discussed above for the temperature), one 
can use bosonization even if some biases are large com¬ 
pared with other characteristic energy scales in the sys¬ 
tem and the problem falls outside the linear-response 
regime. Bosonization thus provides us in many a case 
with a powerful way of addressing strongly nonequilib¬ 
rium transport problems. 

Given the vast array of possible experimental situa¬ 
tions, there are many types of setups to be considered 
and one has to proceed in a class by class basis [6]. We 
shall restrict ourselves to the situation in which the leads 
are Fermi liquids in different equilibria (as is typical of 
Landauer-style setups; cf. Refs. 7 and 8) and the nonequi¬ 
librium situation is confined to a zero-dimensional system 
(i.e., at most a Hilbert space with just a few degrees of 
freedom) adjoining them via point contacts. In particu¬ 
lar, we shall focus on the important example of tunneling 
junctions of different types and steady-state conditions. 


II. CASE STUDY: A SIMPLE JUNCTION OUT 
OF EQUILIBRIUM 

How to set up the bosonization formalism so that it re¬ 
mains valid under nonequilibrium conditions is a delicate 
procedure, often used but not so very often discussed in 
the literature. We shall focus on the important case of 
nonequilibrium steady states and we will expand on the 
details via a particular example that will also serve to 
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highlight certain subtleties whose resolution will be the 
main focus of this paper. 

We choose to study the problem of a point-contact 
junction for noninteracting spin-1/2 fermions formed be¬ 
tween two leads separated by a tunneling barrier. This 
situation is captured via a standard tunneling Hamil¬ 
tonian [9]. Provided the interactions in the leads are 
screened so one can describe them with a noninteract¬ 
ing model, the problem can be reduced to one dimension 
via the use of symmetries [10] or, alternatively, intro¬ 
ducing a lattice regularization and applying a Lanczos- 
Haydock recursion [11] which is valid even in the pres¬ 
ence of disorder. One ends up with a semi-infinite chain 
with the junction or impurity attached to its bound¬ 
ary. In a continuum description (after introducing an 
appropriate bandwidth cutoff and linearizing the spec¬ 
trum) we have two degrees of freedom which can be 
called the incoming and outgoing electrons, that move 
towards or away from the boundary, respectively. These 
are defined for r > 0 and obey the boundary condition 
= 0,t) — (r = 0, t); we omitted internal in¬ 
dexes for brevity. One can introduce the new operator 
■0^ (r, t) = 9 (—r) (—r, t) + 9 (r) 0^^^ (r, t) defined in 

the whole axis, where 9 (r) is the Heaviside step function. 
These are now chiral fermions and the junction problem 
has been mapped to two chiral-fermion leads adjoined 
at a point (see Fig. 1). These standard transformations 
are commonly used in the context of (quantum-)impurity 
problems [12, 13] (see also Refs. 14 for a related descrip¬ 
tion in the context of numerical renormalization). No¬ 
tice that we could equally well exchange 0,^^^ and tplut 
while defining 00 which means that we have the free¬ 
dom to work with either right or left movers and which 
type to choose is a matter of convention. In what fol¬ 
lows, after the mapping to a one-dimensional model, we 
will refer to the spatial axis as the x axis. The open- 
circuit (open-junction) boundary conditions now read 
0’!’ {x = 0~,t) = {x = 0+,t). 

A. Setting of the Problem and Direct Solution 

In Hamiltonian language, the model we shall focus on 
is given by 

’ ( 1 ) 

where the Hamiltonian (densities) describing the leads 
and the Hamiltonian for the tunneling across a barrier 
modeled as a quantum point contact between the two 
leads [17], respectively, are 

Hi = VF'iplgix,t) {-id:,)'ipae{x,t) (2a) 

Htnn = -It i’ii (0, t) (0, t) . (2b) 

Here ipai{x,t) are spin-I/2 (a =t,i) chiral fermions in 
the Heisenberg representation that are obtained after 



FIG. 1. Schematic representation of the mapping procedure. 
When each lead can be described as fermions on a half line 
and there are no interactions (or at least no backscattering), 
one can unfold the space into a full line with only one type of 
movers. These are called chiral fermions and the final setting 
is in many ways similar to that which presents itself naturally 
in setups involving quantum-Hall edge states. By convention, 
one sets the junction at the point x = xq = 0 (with x being 
the vertical axis in the figure). In the present case to be 
studied, the junction is simply a potential barrier modeled 
with a tunneling-overlap matrix element yt. More in general, 
one could have a more complicated tunneling system (such as 
a double barrier with resonant levels in between) in order to 
describe tunneling through nanostructures such as quantum 
dots. The unfolding procedure can also clearly be generalized 
to three-terminal settings and more [15]. 


“unfolding” the two leads as described above. Notice 
that the leads are modeled with an exactly linear dis¬ 
persion and Up is the Fermi velocity. The tunneling ma¬ 
trix element, yt, that characterizes the barrier is taken 
to be energy independent and the local fields at a; = 0 
are understood in the sense of the local-action formal¬ 
ism used for the calculations as explained below. We 
took here yt to be real for notational simplicity (though 
that is not necessary and we will write more general ex¬ 
pressions later). Going beyond the standard (physically 
motivated) regime, we will consider jytj to be arbitrarily 
large. To fully define the physical situation, we still need 
to describe the (nonequilibrium) state of the system. We 
assume now that at a much earlier time the connection 
between the two leads was established and that there is a 
battery keeping a constant chemical-potential difference 
between the two leads. Let us call fj-e the chemical poten¬ 
tial of lead £ = L,R= such that = eV and the 

full potential drop takes place in the junction region (see 
Fig. 2 for a sketch of the physical configuration). The in¬ 
formation about these chemical potentials will enter into 
the distribution functions for each lead. Under these con¬ 
ditions we know, by design, that the system would have 
reached some nonequilibrium steady state [18] and can 
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FIG. 2. Schematic depiction of the setting in which two Fermi 
seas kept at different chemical potentials with a difference 
given by eV = pl — Pr are connected via quantum tunneling 
across the potential barrier that separates them. The tun¬ 
neling across the classically forbidden region is modeled by a 
tight-binding matrix overlap 7 t that can be taken to be energy 
independent in certain cases (in particular, the characteris¬ 
tic energy scale of dependence of the transmission coefficient, 
| 7 t|^, has to be much larger than both eV and k^T [ 8 ]). Here 
we assume that the barrier region does not allow for inter¬ 
nal states; a situation when that happens will be discussed 
elsewhere [16]. 


be described with a time-translationally invariant action 
(this is only important for the particular solutions we 
discuss below, but not for the more general conclusions 
that we reach). 


1. Transport Characteristics 


Green’s function of the junction [20] , 
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(4) 

where 

sn = St (w) = 

1 

- 2 /( 

' UJ-flg \ 

, Tt 1 

= tanh 

^ and Tt 


is the temperature of each lead; / {x) = l/(e’’’ -b 1 ) is the 
Fermi function. The current can be computed according 
to 


A A/^ ? ? 

I = dt— = -[H,AN] = - [Htun, Nn - iVn] 

= (0, t) (0, t) ^ (Ji = 0) (5) 


(we shall follow the convention in which sums over vary¬ 
ing indexes are implicit). Thus, restoring the complex 

conjugate tunneling amplitude 7 ^, we have / = ^/^ = 

— / [Tt^RL ~ 7 t ^LR ]: where the expectation value is 
evaluated via a choice of appropriate Green’s functions in 
the Schwinger-Keldysh formalism. Next one can proceed 
to invert the inverse of the Green’s function and find the 
necessary expressions for the integrand, 


7tGR+ 


7t*G: 


- + 
LR 


2 \t\^ (sl — 5r) 

(i + wf 


( 6 ) 


One of the reasons for defining the problem as we did is 
that it is amenable to an exact solution. How to find the 
transport characteristics is well known and there exist a 
number of standard ways of going about it. So we shall 
be brief but give a complete summarized account in order 
to highlight notations and conventions. Our approach is 
to integrate out the degrees of freedom in the leads that 
are not directly active in the tunneling process (i.e., away 
from X = 0 ) and thus derive a local action for the problem 
[2]. When doing that, one obtains the diagonal matrix 
elements as momentum-space integrals of the two-point 
Green’s functions regularized as principal-value integrals. 
Notice this is consistent with the normal-ordering and the 
regularization of the diagonal terms in the action, which 
are needed for the bosonization treatments that are the 
focus of this work. One can capture the non-equilibrium 
situation by using a standard Schwinger-Keldysh formal¬ 
ism (see Refs. 19 for examples). At the moment we ne¬ 
glect the spin, which will just give a factor of 2 at the 
end, and choose the following spinor basis in Keldysh 
space (here we follow the same notation as in Ref. 20, 
but we reorder the basis): 

^ (3) 

where the index k labels the Keldysh-contour branch fol¬ 
lowing the ‘minus-means-forward’ convention [21]. Let us 
define "ft = 2vFt and use the result for the local inverse 


More explicitly, and including now the spin-degeneracy 
factor, we have 


I = 


4|tr 


4 If 

(l + lf)' 


J [SL (w) - Sr (w)] ^ 


(7) 




[/r (w) - /l (w)] duj. 


This expression has a standard form and is intuitively 
appealing, as the integrand selects a window (smeared 
by the temperature) that is 1 in the frequency interval 
between the two chemical potentials and zero outside of 
it. The integral can be carried out in general, but we 
will be taking the zero-temperature limit (Tl = 7r = 
T —i 0) for simplicity. In that limit, both become sign 
functions (/^ become step functions) and the integration 
is trivial: 


I —^ 

T-fO 


4|f eK 

7r(l-k|t|^) 


( 8 ) 


This gives the particle current, and as always one needs 
to multiply by (—e) to get the electric current instead. 
The result is standard [17] and the fact that the response 
is exactly linear in V to all orders is a property of the 
linear spectrum of the model. 
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B. Bosonizing in the Steady State 

On the one hand, to study a problem using bosoniza- 
tion, one of the first things to do is to factor out the 
fast modes [1]. On the other hand, to study a problem 
in which a finite voltage bias is present, one of the first 
things to do is to introduce it into the calculations (for 
instance, via a careful treatment of the interaction pic¬ 
ture [22]). Here, we need to take care of both things, so it 
is better to discuss them in the more formal unified lan¬ 
guage of gauge transformations. (See, though. Ref. 23 for 
an approach in terms of scattering states that provides 
an alternative to ours but is ultimately equivalent [24].) 

As discussed already in the introduction, bosoniza- 
tion is a rewriting of the excitation spectrum in terms of 
bosonic degrees of freedom. As such, it does not capture 
the information about the reference state or ground state 
(which is a Fermi-Dirac sea of noninteracting fermions). 
Technically, one would say that what one knows how 
to bosonize is the normal-ordered Hamiltonian in which 
the vev (vacuum expectation value) has been subtracted. 
The type of nonequilibrium situation we are considering 
here presents thus a problem, because one knows in prin¬ 
ciple how to normal-order for each lead, but only in an 
open-junction configuration. The subsequent inclusion 
of the tunneling term constitutes a delicate task. A sys¬ 
tematic way of carrying this out starts by using time- 
dependent gauge transformations to map the finite-bias 
problem into a zero-bias one but with explicitly time- 
dependent couplings. 

Let us first switch to Lagrangian language (to fully 
capture the effects of a time-dependent gauge transfor¬ 
mation), in which the system is described by the La- 
grangians (densities): 

t) 'ipai {x, t) - 

= (*^t + WF^o;) V’rrt (x, t) , (9a) 

Atun = --fftun = 7tV'L (0> 0 P’al (0. • (9b) 

We can now make the following field transformation 
tpcrtixjt) = (notice that some authors 

follow as an alternative a prescription of including a time 
dependence related to the lead chemical potentials into 
the respective Klein factors when bosonizing the model 

[25] ). This gives 

(x, t) {idt + Ppat {x, t) , (10a) 

Ltun = (0^ t) (0, t) 

= ( 0 : t) ( 0 . t) + 

-k (0,0'0 <tL (0, <) (10b) 

(where, in the last line, we restored explicitly the com¬ 
plex conjugate yC). An important point is that now the 
distribution functions in the Keldysh action do not con¬ 
tain information about the chemical potentials any longer 

[26] . Next we subtract the vev’s of each lead which, by 


assumption, are the same as those in the absence of the 
tunneling term (this is where the Landauer prescription 
[8] enters the calculation) and drop the now ineffectual 
chemical potential terms. For a noninteracting problem 
this is equivalent to factoring out the fast oscillations 
in each lead according to an additional field transforma¬ 
tion: V'trf [x, t) = {x, t), with kp = fii/vF for this 

linear-dispersion case, and then subtracting the same (in¬ 
finite) constant for all leads. So we are naturally left with 
the normal-ordered formulation of the problem, 

= ■ '^ii {x, t) {idt + ivFdx) 'ipai {x, t) : , (11a) 

Ltun = (0, t) PpaF (0, t) + 

-Fe"*®'^‘7t^iR(0>0^<TL(0,<) . (11b) 

At this point we lost all the information about any 
absolute-energy reference, but we still have the infor¬ 
mation about the potential drop encoded in the time- 
dependent phase of the tunneling term (cf. Fig. 2). Given 
the infinite-bandwidth setting, we are also in a situation 
in which space is naturally to be regarded as half filled. 
Now one is ready to bosonize the problem following the 
standard procedure. 

1. Abelian Bosonization and Standard Transformations 

The Abelian-bosonization recipe is by now textbook 
material [1, 2, 5, 27] and there is no need to present 
the details here. There exist though a number of differ¬ 
ent conventions, which can bring in some confusion at 
times. Our notation and conventions follow closely the 
review article in Ref. 4 (which in turn is based on the con¬ 
structive presentation given earlier by Haldane [3]), with 
the only difference of factors of that are needed 

in order to have a more standard normalization for the 
real-space Fermi-held anticommutators [28]. 

In order to bosonize we go back to the Hamiltonian 
formulation 

(^- 0 i-'ivFdx) {x, t) : , (12a) 

-fftun = (0> ^<tR (0, t) - 

- e"“^‘7t ^Ir (0, t) Pp^F (0, t) , (12b) 

and we proceed to bosonize according to 77°, which is 
akin to working in the interaction picture (with TLtun 
taken as the interaction term) [29]. We shall follow the 
bosonization prescription [30] 

Pp^t {x, t) = {t) , (13) 

V27ra 

where the Fg-i (t) are the so-called Klein factors and a 
is a short-distance regulator [3]. We shall not include 
subleading 1/L corrections in the bosonization formulas, 
because inhnite size is the appropriate limit for a descrip¬ 
tion of the leads in a Landauer-style transport setup to 
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describe a steady state; as a bonus, this keeps formulas 
shorter. In terms of these bosons the Hamiltonian den¬ 
sity for the leads can be shown to take the usual form, 

= = ^ . (14) 

E cr—L,R 

One of the main advantages of the bosonic descrip¬ 
tion is that with it one can more easily recombine de¬ 
grees of freedom in order to, for instance, separate the 
effects of charge and spin dynamics (phenomena such 
as spin-charge separation are thus very naturally de¬ 
scribed with the use of bosonization). Using the stan¬ 
dard, physically motivated, rotated boson basis (pai = 
\ {<t>c + o'(j)s + + cricfisi), where a,i = ±1 when enter¬ 

ing as multiplying factors, the noninteracting Hamilto¬ 
nian density retains its quadratic form, 

^ : (15) 

iy—c,sj,sl 


we will need the following ones: 



(17e) 

FlRFiL=F!Fsi, 

(17f) 

F^rF-^r = FiFsi , 

(17g) 

fIrF^r = F^Fi , 

(17h) 


(where the last two are simply the Hermitian conju¬ 
gate of the first two). Notice that, as one should have 
expected by looking at the boson fields and compar¬ 
ing Eqs. (16a) and (16b), the right-hand sides involve 
only the lead and spin-lead Klein factors. The tunneling 
Hamiltonian density can thus be further rewritten refer¬ 
ring only to the ‘physical’ sectors. One can then undo the 
steps of the bosonization procedure and debosonize (also 
called reverse bosonization or refermionization) in order 
to arrive again at a problem written in terms of Fermi 
fields. Using the standard debosonization prescription, 
■01/ {x,t) = (t) which parallels the one 

we used for bosonizing in the first place, we arrive at 


and, as usual, the Klein factors drop out from these 
terms. We shall refer to these ‘physical’ sectors as charge, 
spin, lead (or flavor), and spin-lead (or spin-flavor), re¬ 
spectively. We will see how they naturally reorganize the 
information about the physics of tunneling transport. 

Let us now bosonize the tunneling term, rotate the 
bosons into the physical sectors, and make some standard 
simplifications (the sum over a is implicit and the fields 
are evaluated at a; = 0): 


Htur, — —t 


eVt 7t p,t p 

_ X r ra-RC 




27ra 
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2'Ka 


O' 


_peVtJi^p^ p 

27ra 


7t p-j p^_^pi4'l+cr<pBt) 
2'Ka 


(16a) 


(16b) 


To proceed further we need to take care of the mapping 
of Klein factors. We anticipate no subtleties coming from 
these, but we carry out a careful treatment nevertheless 
so as to show that explicitly. The most rigorous way 
to proceed is by identifying relations between different 
bilinears of old and new Klein factors, and fixing the 
four arbitrary phases that appear [31-33]: 


= fIf! , 

(17a) 

II 

(17b) 

F^rF^l = fIf^ , 

(17c) 

II 

(17d) 


njO 

tty 

Htun 


--: {x, t) {-ivpd^c) '0R {x, t) : , (18a) 

-^‘7t^z(0,t) + e-“^‘7t>J(0,i)' X 

X l-ipsi (0,t) - (0,f)j . (18b) 


We find that the tunneling term involves only the lead 
and spin-lead sectors, while the charge and spin sectors 
have decoupled from the tunneling process. 

The new problem, defined by = f 91° -h Htun, can 
now be regarded as arising from an original problem with 
the voltage acting as a chemical-potential shift of the lead 
fermions only {v = 1). In other words, if we consider the 
problem given by 


91° == ipl (x, t) {-ivRdx) V’r {x, t) 


[btV’i (0, t) + 7 *tpj (0, t) (0, t) 


Htun — — 


(19a) 


(19b) 


where the chemical potential is set as = — (eU) and 
is zero for all sectors v ^ I, this can be connected with 
the debosonized problem of interest following equivalent 
steps to those we presented above via the combined trans¬ 
formation {x,t) = ix,t). Moreover, this 

‘parent’ problem can be seen to be unique (i.e., there is 
only one way to eliminate the time dependence from the 
tunneling term by reintroducing chemical potentials into 
the problem). 


C. Indirect Solution using Conventional 
Bosonization-Debosonization 


All the rest of the Klein-factor bilinear relations can be One of the goals of a bosonization-debosonization pro- 
derived from these. In particular, in order to simplify gram (BdB for short), as exemplified above, is to achieve 
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a simplification of the problem at hand that would not 
be so easy otherwise. (There could be other alternative 
or additional motivations for bosonizing, such as carrying 
out a renormalization-group analysis that is more easily 
done in the bosonic language; see Ref. 1 for examples.) 
Indeed, transformations like the one introduced by the 
simple rotation of the bosonic basis would be hardly ev¬ 
ident if one were to express them directly in terms of 
the old and new fermions instead. The example that we 
picked is special, because we are able to solve it exactly 
already in the original formulation and even in an out- 
of-equilibrium setting. However, the BdB program is, in 
most other cases, crucial for simplifying the problems and 
being able to hnd solutions either exact or approximate. 

In the case of our simple junction problem, the BdB 
program does indeed show some apparent simplihcations. 
A simple glance at the final form of i?tun shows that 
only the lead and spin-lead sectors are involved in the 
transport while the other two sectors {spin and charge) 
do not participate. This provides a certain economy of 
description that we will discuss further below. For now, 
our immediate goal in this section is to recompute the 
I-V characteristics of the junction. 


1. Re-calculation of Transport after Conventional BdB 
We need again the operator expression of the current, 


but now in terms of the new fermionic degrees of free¬ 
dom. One can translate it from the expression we gave 
above [see Eq. (5)] using BdB or, equivalently, it can be 
recomputed directly in terms of the new fields: 


fa*" ■ 
7 = a—= ■ 


H, 


AN' 




tun, Nif—i 


( 20 ) 




V’L ( 0 ,t) - ' 0 s/ ( 0 ,t) 7 t 0 / ( 0 ,t) - ( 0 ,t) 


Notice that this time the spin degeneracy is already in¬ 
cluded implicitly in the formalism. Thus, / = ^/^ is 
given as 

I = -*7t (( 0 I/ (0, t) ijji (0, t)'j - {ipsi (0, t) 0/ (0, t))^ -b 

+ * 7 t ((0I/ (0,00/ (0,0) - ( 0 s/ (0,00/ (0,0)) 


Next we calculate the necessary Green’s function ele¬ 
ments using the same procedure as in Sec. IIA I. How¬ 
ever, this time we need to introduce a Nambu structure 
due to the presence of anomalous processes in i7tun- As a 
result, we adopt the following spinor basis (including also 
the Keldysh indexes and with the frequencies restricted 
to the positive semiaxis only in order to avoid double 
counting): 


«'(a;) = (0; (w) 0;+(w) 0]' (w) 0/+(a;) 0,,; (w) 0jj (w) 0^, (w) 0^+(w))^ . (21) 

We write the local inverse Green’s function of using the fact that all non-equilibrium Green’s functions (i.e., 
advanced, retarded, and Keldysh components) are diagonal in the Nambu basis. The only change required for the 
time-reversed Nambu component, as compared with the time-forward one, is to dehne Si, = tanh ^ given 

in the argument of the spinor (and we will be taking the temperature to be uniform, Ti^ = Temp)- Including also the 
contribution of i?tun, the local inverse Green’s function for the junction is thus given by 
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We invert the matrix, identify the relevant matrix ele¬ 
ments, and replace them into the expression for the cur¬ 
rent. After some algebra one gets 


/ = 


l^f 

(l + |tf) 



[si (w) 


Si (w)] 


duj 
^ ’ 


(23) 


where t' = 2t. The integral can be done in general, but 


in the zero-temperature limit reduces to 



We see that the result we obtained for the current 
shows several discrepancies from the one in Sec. HA I. 
Such differences need to be understood. 
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III. THE NONEQUILIBRIUM TRANSPORT 
PUZZLE 

We have carefully chosen the nonequilibrium junc¬ 
tion problem so that it meets all the requirements for 
bosonization to be an exact operator correspondence be¬ 
tween fermions and bosons (cf. Ref. 4). All the transfor¬ 
mations we carried out are thus rigorous and the discrep¬ 
ancy between the results of Secs. IIA 1 and IIC I is not 
only unexpected but also unwelcome. There has to be 
an inconsistency somewhere and, given that the result of 
the direct solution is standard and can be reobtained in 
a number of alternative ways, everything seems to indi¬ 
cate that the problem has to be with the indirect solution. 
Moreover, the actual transport calculation of the indirect 
solution proceeded in a very similar way to the case of 
the direct one. As a result, the reason for the discrepan¬ 
cies is likely not in there, but in the preceding BdB-based 
mapping used to rewrite the junction problem in terms 
of the new fermionic degrees of freedom. 

Before furthering the analysis, let us first catalog the 
discrepancies between the two solutions: 

1. To match the solutions one needs to arbitrarily cor¬ 
rect the tunneling matrix element of the indirect 
solution by a factor of 2 (namely, t' i—?> t) in order 
to make it look closer to the exact direct solution. 

2. There is a overall factor of 4 difference between the 
two solutions (the indirect solution would need to 
be multiplied by 4 to match with the direct solu¬ 
tion). 

3. There is also an additional factor of (1 -I- |t|^) in 
the numerator of the indirect solution that cancels 
one power from the denominator and introduces a 
further discrepancy with the exact direct solution. 


These three discrepancies are present no matter which 
method we use for the final transport calculation (they 
all yield the same result). We highlighted them by look¬ 
ing at the zero-temperature limit, but it is easy to see 
that they are also exactly the same at finite temperature. 
Additionally, very similar discrepancies can be seen to be 
present in equilibrium thermodynamic calculations using 
a Matsubara formalism (see the Appendix). Thus, the 
puzzle is not restricted only to transport, but it is more 
evident in transport calculations. 

Motivated specially by the third entry from the list of 
discrepancies, one could imagine expanding the results 
of the direct and indirect solutions in powers of t. It is 
clear that big differences will show up as soon as one goes 
beyond leading order in the tunneling matrix element for 
both calculations. We therefore expect to be able to gain 
some insight by studying the problem using perturbation 
theory in t. 

A. A diagrammatic diagnosis 

Let us start by setting up a dictionary for processes 
allowed by the different vertexes in Htun- There are four 
of those, given by the two possible spin orientations and 
the two possible directions of tunneling. Since our BdB 
program rests neither on the SU (2) invariance nor on 
the hermiticity of iLtun, we can, in principle, set the four 
corresponding matrix elements to different constants and 
thus individually trace each process thorough the BdB 
procedure to construct the dictionary given in the table 
below. Alternatively, one can construct the dictionary 
by looking at the changes operated by the different graph 
vertexes on the fermion numbers of the different sectors 
(which is essentially the construction that is used to iden¬ 
tify the different Klein-factor bilinears [31]). The trans¬ 
lation between the fermionic structure of the vertexes in 
terms of ‘old’ (original) and ‘new’ fermions is thus given 

by 


Simple-junction Graph-vertex Dictionary 

Original Fermions I j New Fermions 


^IlV’tr t R 
V’lLV’tRjiR 


tR 

;r 

tL 


i’li’si 


si 


si 

si 


I 

l 


Notice that the second two lines are the Hermitian conjugate of the first two. We can refer to them as (i)-(iv) from 
top to bottom. Now in order to calculate the current we need to find the fully dressed vertexes. We can proceed to 
dress them by carrying out a perturbative expansion in LItun (the Keldysh structure is not important for the present 
argument and will be suppressed for the sake of clarity). 











Let us consider, for instance, the dressing of vertex (i) in terms of the original fermions. It proceeds by alternating 
vertexes (i) and (iii) at different orders of expansion. Up to third (the first nontrivial) order we have 

t L —t R = t L ► • ► t R + 

+ tL - ^ -•-►-•- ^ - tR- + 

, tR tL 


which in terms of the new fermions translates according to our dictionary into 

I ^ —►— si = I ^—• ► si + 

+ I - < -•-►-•- M^ - si + 

, si I 


For spin-down, the diagrams in terms of the original fermions are exactly the same with the obvious label replacement 
(t) (i)- This corresponds to alternating vertexes (ii) and (iv) at different orders of expansion. After translation 

to the new-fermions language one just changes the labels according to (j, L) —)• {sl) and (j, R) —>■ {1), but this time 
the arrows of the fermion propagators stay unchanged (no anomalous processes are involved in this case, exactly the 
opposite from the example above with spin-up). 

Difficulties arise when we start directly from the new-fermions language and proceed to dress the vertex in question. 
This is so because we have additional (and, we shall claim, unphysical) ways of introducing contractions. Consider, 
for instance, again the case of vertex (i). One would proceed to dress it as follows: 


I ^ —►— si 




-f I 
-f I 
-f I 

-f I 

+ .. 



si + 
si + 
si -f 
sl + 


The four third-order processes correspond to vertex insertions (i-iii-i), (ii-iv-i), (i-ii-iv) and (ii-i-iv), respectively; which, 
according to our dictionary, translated back in terms of the original fermions read as follows: 

t L —t R = t L ► • ► t R + 

-btL -►-•-►-•-►-•-►- tR + 

tR tL 

-|- t L -•-► t R T 

tR tL I tL 

-btL -»-•- > -•-^»- iL + 

tR I tL tR 

-I- t L -►-•-►-•-►-•-►- t L -b 

, tR|tLtR|tR 


The last three contractions are not allowed in the 
original-fermions framework, as they require spin flip and 
some even L R exchange (as indicated by the inner 
labelings). Moreover, they do not even dress the cor¬ 
rect vertex (as indicated by the outer labelings). From 
a practical point of view, one may notice that while we 
deal with four distinct types of original fermions (f L, 
I L, t R and j, R), we deal with only two types of new 
fermions (/ and si ). We conclude that the more compact 


description achieved after the BdB-based mapping intro¬ 
duces the possibility of spurious processes that should 
not have been there. These are processes that mix ver¬ 
texes (i) and (iii) with vertexes (ii) and (iv), which in 
terms of the original fermions cannot happen due to spin 
conservation. This clearly hints at the possibility that, in 
the new-fermions framework, the spin sector should not 
really be decoupled after all. 
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IV. CONSISTENT APPROACH TO 
BOSONIZATION-DEBOSONIZATION 

We need to revisit the transformations in the BdB- 
based mapping used above, with the goal of finding the 
source of the discrepancies with respect to the direct cal¬ 
culations. In particular, one needs to be careful about 
the fact that the tunneling term is not normal ordered 
(since the procedure of subtracting the vev is not well de¬ 
fined for the processes in i?tun for they are not diagonal 
in fermion ‘internal indexes’). 


A. Keys to Consistency 


We proceed to study again the bosonization of the tun¬ 
neling term but taking care of not combining exponen¬ 
tials. If we start from Eq. (16a) and perform the change 
of basis for the bosons, we arrive at 


ff _ „ieVt 7t pt p i{<pa+cr 4),-<f)I -cr 4)si)/2 

xitun — e c tCctRc a 

ZTTO 


X e 


-i{4>a-i-cr4>s-\-4>l+<^/>Bl)/'2_ 


„-ieVt 7t t7.t P i(4>a+(y4>s+4>l+’^<t>sl)/'2 


X e 


-i{4>a+(74>s-4>i-'y4>si)/‘i 


These new objects (to be defined and discussed more in 
detail below) can be interpreted as particle and hole den¬ 
sities for new fermionic degrees of freedom in the eharge 
and spin sectors. They shall be considered in their ‘eigen- 
basis’ and they have eigenvalues 0 or 1 and 1 or 0, respec¬ 
tively and correspondingly. This is the central result of 
the consistent way to debosonize and the (almost) final 
form of the tunneling Hamiltonian is 

Htun = 

It can be easily seen that the inclusion of the n factors 
naturally avoids the mixing of graph-vertexes (i) and (iii) 
with graph-vertexes (ii) and (iv), exactly as was con¬ 
cluded to be necessary in the diagrammatic discussion 
of the previous section (Sec. Ill A). Notice also that, in 
the same vein, these factors also stop us from being able 
to rewrite Htun in terms of Major ana-fermion combina¬ 
tions. 

In the next two subsections we provide some additional 
rationale, but those readers that want to skip some of the 
technical discussion can jump ahead to the last subsec¬ 
tion of this section (Sec. IV D) and see how we are now 
able to recover exactly the results of the direct calcula¬ 
tion (which can be taken as a pragmatic justification for 
the procedure). 


We will now, on the one hand, combine the exponentials 
in which the bosons appear with the same sign (we are 
prompted to do this by a study of the corresponding op¬ 
erator product expansions, OPEs, and by the consistency 
with the mapping of the Klein factors [31, 32]). On the 
other hand, we will be careful not to combine the ex¬ 
ponentials in which the bosonic exponents appear with 
opposite signs (prompted by the suspicion, from our per¬ 
turbative analysis, that the v = c,s sectors should not 
completely decouple from the tunneling process). We 
will discuss the charge and spin sectors carefully momen¬ 
tarily; for now we debosonize in the lead and spin-lead 
sectors only (using the same prescription that was intro¬ 
duced above). The tunneling term takes the following 
form (all the fields are evaluated at a; = 0 and at time t): 

Htun = 

which is the same as before but with the addition of the 
extra exponential factors. 

A pragmatic way to proceed in order to debosonize in 
the charge and spin sectors as well is by replacing the ver¬ 
tex products by lattice-like fermionic densities according 
to the prescription 

^ ( 25 ) 


B. Matters of Regularization 

The exponentials of bosonic fields of the type 
are central objects in the bosonization formalism known 
as vertex operators. The bosonization prescription tells 
us that ipl oc e'^^" (with A = 1) while normal-ordered 
densities are bosonized according to : 

The consistency between these two prescriptions can 
be checked by bosonizing the non-normal-ordered case, 
V’JV'i/ = e®'^‘'e“®'^‘"/27ra, and expanding the right-hand 
side by using known results for the OPEs of vertex oper¬ 
ators [4]. 

However, when bosonizing, oftentimes our aim is to 
change basis from the spin-&-lead (or spin-&-flavor) 
states to a basis that separates physical sectors {charge, 
spin, lead and spin-lead) because some of the physics will 
simplify by doing that (this is the transformation that 
we performed in Htun)- Proceeding formally for each 
vertex operator of a density operator [using (x) = 
where the label the physical sectors and 
we absorbed minus signs that are not important for this 
part of the discussion], we have 

(26a) 

n 

[1 + a d(j)cri + ...] ~ 1 -f ^ ^ Q /2 + ... (26b) 

n 

1 + Sufyi -h ... ~ 1 -f ^ ^ 6ni,^l2 + ... 


(26c) 
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where in the third line we introduced lattice-like density 
fluctuations, Sriai = ad(f)a-e and = ad(j)^^^ to stress 
that they need to be small in order to connect to the first 
line. (A standard view is to treat a as a control parameter 
for the expansions in the second line.) Therefore, these 
transformations are consistent if bosonization is treated 
as an expansion around a half-filled ground state (in a 
real-space picture). While the bosonization identities are 
precise, some manipulations might not hold when the 
deviations from the local half-filled state are large. If 
a particular problem, as is the case of some transport 
problems like the one that we are studying, forces us to 
consider large 5n fluctuations, then we need to proceed 
with caution while expanding. 

One solution is to expand around a different state, 
which can be achieved via a linear transformation. Con¬ 
sider the following vertex OPE at some x = xq (the po¬ 
sition of the junction or impurity) and treat a as an ex¬ 
pansion parameter (not necessarily small) [34] : 

« 1 + a + ... = 1 + (l + a + ... 


as will be proven below by working to all orders without 
resorting to OPEs (see Eq. 33). 

There is a delicate point regarding the proper nor¬ 
malization (or scaling of the coupling constants) of non- 
normal-ordered terms as those in i?tun- This is more 
easily understood considering vertex products diagonal 
in internal indexes. Using the consistent identities de¬ 
rived above, we can proceed as follows: 



(30a) 

- y/16 

(30b) 

y n 

(30c) 

n 

(30d) 


(30e) 


n 


This serves as a definition of a shifted set of bosons, 
(pa-i, which are used to expand around a differently 
filled state (unit-filling in this case) and need to obey 
dpae = dpa-e — 1/a. Reintroducing the x dependence 
from the OPE before taking the a —^ 0 limit [i.e., replac¬ 
ing 1/a ttS {x — xo)], and integrating this relation one 
gets (pa-e {x) = (pai (x) — f sgn (x — Xg), up to an additive 
constant. The new bosons have identical commutation 
relations and OPEs except at x = xg due to the presence 
of these solitonic shifts. 

Expanding around adp^e = 1 is equivalent to expand¬ 
ing around adpo-e = 0 and we can use small-variable ex¬ 
pansions in terms of the latter. For the kind of vertex 
products we are considering (at x = Xg) we have 

« 2 (l -H "^dp^e -h ...) -H ... (27a) 

« 2^1 + adp^e + ... (27b) 

~ ; (27c) 

where to get to the second line we used a Taylor expan¬ 
sion for the square root of a binomial (the first two lines 
are strictly equivalent to the order that is given explic¬ 
itly; their connection can be regarded as a sort of partial 
re-summation that is also consistent with a further study 
of other vertex OPEs that we carried out as well). Al¬ 
ternatively, applying the same vertex-vertex OPE, but in 
reverse, to the parentheses in the first line of the equation 
above we have 

This implies that, generically. 


(29) 


where (i) in the first line we shifted the bosons away from 
half filling; (ii) from the first to the second line we did 
a change of basis; (iii) going to the third line we used 
Eq. (29); (iv) from the third to the fourth line we shifted 
the bosons back to half filling and we also distributed 
the overall square root; and (v) finally we redistributed 
the square root between the two vertex operators again 
using Eq. (29). Notice the introduction of the primes (in 
PP) to distinguish this case when the change of basis is 
done with the p's, from the case when it was done directly 
with the original p's. The primes will be dropped when 
a comparison is not being done and the case in point is 
clear from the context (this notational variation is used 
in this subsection only). 

What we found is that if the change of bosonic ba¬ 
sis is done in terms of shifted bosons, then a prefactor 
of 1/2 appears for proper normalization (and we shall 
make this conclusion extensive to nondiagonal products 
as well). This kind of normalization changes, or rescaling 
of couplings, is common in bosonization treatments and 
can often be traced to subtle differences in regularization 
schemes. In particular, shifting the bosons is equivalent 
to acting with so-called boundary-condition changing op¬ 
erators [2, 35], which are a known source for “coupling- 
constant redefinitions” (for another example, also involv¬ 
ing a relative factor of 2, the reader can look at Sec. 2 
of Appendix A in Ref. 31). To summarize our result, 
we should contrast the differences between Eq. (30e) and 
the one we presented at the start of this subsection in 
Eq. (26a). To develop some intuition, let us introduce 
the lattice-like notation 2nv « 1 -|- 5n,j, so that, near 
half filling, the left-hand side is close to 1 and near maxi¬ 
mum filling it is close to 2. We will also use the notation 
•\/2hi, « y/2ji// « vT+~5^ (these will be made more 
precise in the next subsection). The two BdB-mapping 


g*,^/2g-i0/2 _ ; 











11 


relations can then be rewritten as 

2ncre « 'J2ny^ when 2nai ~ 1, (31a) 

n 

~ 2 n V2h,y^ when 2nae ~ 2 (or 0); (31b) 

n 

where we highlighted that they are useful in different 
regimes. Which one, or when each of the two, should 
be used needs to be judged depending on the problem 
that is being solved (and that is part of what we mean 
by a consistent use of BdB-based transformations). We 
argue that the junction problem requires the use of the 
second one, because the physics of tunneling calls for the 
consideration of unit-size particle-number fluctuations at 
the junction {uat = 00 1). 

It is instructive to see how these two different regimes 
(i.e., half filling versus maximum/minimum filling) are 
connected in our formalism to a change of boundary con¬ 
ditions for the new fermions after the BdB-based map¬ 
ping. We start from the continuum boundary condi¬ 
tions, V'L(0~) = V'L (O”*^) for ^11 rrf. After changing 
basis in the intermediate bosonic language of the (/)’s, 
we get to the new fields with (0“) = (O'*") for 

all = c,s,l^sl, as naturally expected. If we do 
the change of basis with the (f)’s instead, the result¬ 
ing boundary conditions are different. In the charge 
sector, from the definition of the c^’s it follows that 
(t>c = 4>c — 7rsgn(a;), (where we went back to xq = 0), 
and thus (j)'^ = (pc + |sgn(a;) = — |sgn(a:). From 

there it follows that ip\ (0“) = —ip\ (0+). For the other 
sectors (i>n ^ c), we simply have (p^^ = <p^^, and thus 
(pP^ = (pv^ + I sgn(x) = (/)y„ -I- |sgn(a;); so that, in a 
different way, we still get that tpl^ (0“) = — ipl^ (0^). 

Remarkably, these antiperiodic boundary conditions 
parallel what Affleck calls “strong-coupling boundary 
conditions” in the context of the boundary-conformal- 
field-theory approach to quantum-impurity problems [see 
Eq. (1.29) of Ref. 13]. The name is because these are 
the type of boundary conditions needed in the strong¬ 
coupling limit of those problems. What these boundary 
conditions actually do is to decouple the band-fermion 
degrees of freedom at x = xq from the rest of the bulk; 
that way they are not tied to half filling (or other) con¬ 
ditions and they are available to couple them (strongly) 
to the impurity. In our case, we shall in general need 
those degrees of freedom to be available (even if there is 
no impurity) to participate unrestrainedly in transport 
situations. 

We shall thus refer in our context more generically to 
consistent boundary conditions (CBCs). These depend 
on the problem at hand and in the particular example 
studied here they turn out to be antiperiodic boundary 
conditions. Notice that the need for a factor of 1/2 as dis¬ 
cussed above can be seen as the practical manifestation of 
the boundary conditions that were (implicitly) adopted. 
Let us also mention that the use of CBCs does not mod¬ 
ify the form of the kinetic part of the action (i.e. when 


rewriting Eq. 15 in terms of </ —>■ (p'). The solitons that 
we introduce with (p' will induce in 'HP additional slips of 
27r localized to a length scale of a around Xq, but since 
that is the limit of length resolution and the bosonic fields 
are compact with radius 27r, those contributions consis¬ 
tently drop out. 


C. Tunneling of New Fermions 

It is now a matter of a delicate but ultimately simple 
replacement to finish the debosonization of iLtun in the 
charge and spin sectors. We (re)introduce the following 
definitions (all fields are at ccq = 0 and time t): 

V2h+, = (32a) 

= + a d(pc,s = 

2'Ka 'ipl,s'ipc,s = , 

V2fi- = e-*‘^=/2e*-^=/2 (32b) 

= — ad(pc,s = 

2'Ka ipsp’l = \/ 2n7 . 

Thus fi+g (hj) are simply the square roots of the par¬ 
ticle (hole) density of charge or spin fermions ‘at the 
site’ of the junction. For a physical picture, one could 
think of them as corresponding to a single lattice site 
after a lattice discretization with tto as the lattice con¬ 
stant; even though that is not the type of regulariza¬ 
tion adopted when bosonizing [3]. One can explicitly 
check consistency by calculating their squares via the 
equal-time, full operator product. For that we need to 
write the bosons in terms of their creation and annihi¬ 
lation components, (p (x) = (x) + (p (x), which obey 

[p{x),p^ (x')] =-In (l-e-T K—')+“]); (see Ref. 4 

for the notational convention to point-split the product 
and normal-order the vertex operators): 

givHx) ^^ip(x) [x') ^ 

X QpvHx)^rip(x)i^x') iip[x') 

_ l-e~^K^-^')+°] \ 

\J 2L y 1 _ e“T[*0'-a:)-|-a] j ^ 

X \/l — i^') 

L 

— lg*b(a;)g-jb(a;') 

2 

= n, 





( 33 ) 
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where (j) stands for either ±(/)c or ±(j}s- 

The squares of the h’s have the properties that 
{nf g) = nf g and = 0; we shall refer to these as 

idempotence and co-nilpotence, respectively (notice that 
if an operator on a finite Hilbert space is idempotent, 
one of its square roots is the operator itself). In ad¬ 
dition, the sum of their squares resolves the identity, 
+ ^cs = 1- Thus, as we will see below, they can 
be consistently assigned the eigen-expectation-values 0 

or 1, as if ^(n^. 

Let us introduce the notation yto- = 7 th+hJ = 
7 t /2 (where the factor of 1/2 at the end 

is included for a coupling-constant rescaling in accor¬ 
dance with our discussion in the previous subsection). 
The consistent form of the tunneling term is then more 
compactly rewritten as 


one (as bilinears). As a result, any connected perturba¬ 
tive expansion in Htun does not involve the charge and 
spin sectors and the can be treated as c-numbers, 
(restoring the Gaussianity of the problem). Due to global 
gauge invariance for each lead, the final expressions in¬ 
volve always the squares of the h’s and can thus be simpli¬ 
fied thanks to their idempotence and co-nilpotence. An 
alternative equivalent calculational procedure is to set 
the h’s to their different eigen-expectation-values, to do 
the calculation, and to trace over all such values (not av¬ 
erage over, because they are not exactly conserved quan¬ 
tities). This second path is shorter and makes more ex¬ 
plicit the connection with the direct solution. 

The result one gets for the I-V characteristics, by fol¬ 
lowing the procedure outlined above, is what one was 
hoping for: 


Htnn = 

- e—• (34) 

And we can finally gauge out the applied voltage from 
the explicit time dependence, as we discussed already 
for the conventional procedure, by using ipi (x, t) = 
This gives 


= - 


httV’i - huV’/ btiV’i - 


Notice that we are not able to combine the fields into Ma- 
jorana components, as we did in the conventional frame¬ 
work, due to the spin dependence acquired by yto-. We 
see how this time the spin plays a role and starts to show 
up clearly, as expected from our diagrammatic analysis 
of the problem. 



When comparing with the result of the direct calcula¬ 
tion, given in Eqs. (7) and (8), the matching is now ex¬ 
act and all the discrepancies are gone. Namely, (i) the 
t I—> t/2 correction is not required as it happened natu¬ 
rally courtesy of the CBCs; (ii) the spin degeneracy arises 
automatically and the correct overall prefactor arises also 
naturally; (iii) the extra factor of (1 -I- |t|^) in the numer¬ 
ator is not present. 


V. CONCLUSION AND PROSPECTS 


D. A resolution of the puzzle 

We are now ready to recompute the indirect solution 
to the transport problem after debosonizing consistently. 
Revising the expression for the current we find 

I = -ijtt i'lpi (0, t) 'tpsi (0, t)) - ijti (^'iph (0, t) ipi (0, t)'j + 

+ *7tt {^h (0.0 i’j (0, f)) + hti (W (0, t) (0, t)'^ ■ 

We adopt the same conventions as before for the defini¬ 
tion of the Keldysh-Nambu spinor basis and make also 
the same redefinitions of the couplings to factor out the 
Fermi velocity. The expression for the inverse Green’s 
function is like in Eq. (22) with the addition of the 
spin index into the tunneling terms (which is straightfor¬ 
ward, since all the Nambu-off-diagonal components ac¬ 
quire cr =t while the Nambu-diagonal components go 
with a = l). 

It should be remarked that unpaired fields do not 
enter in the tunneling term and appear only in the kinetic 


By focusing on a case study in which bosonization 
is rigorously applicable and, not less importantly, ex¬ 
act calculations are possible and enable detailed com¬ 
parisons, we were able to uncover some subtleties of the 
bosonization-debosonization procedure that had quite 
strong implications. Besides directly comparing the 
mathematical expressions as we have been doing, it is 
instructive to compare the two results graphically. To 
that end, we plot in Fig. 3 the two indirect solutions 
for the differential conductance (G = dl/dV) computed 
conventionally and consistently. As expected, the two 
results only agree in the limit of t —> 0. Expanding for 
small t, the two results coincide to order O (t^) and start 
to disagree in the coefficient of the term (with the con¬ 
ventional result being larger by a factor of 2). This is as 
expected from the diagrammatic analysis. Let us remark 
that, for t > 1, the differential conductance computed 
conventionally not only lacks the t ^ 1/t duality of the 
exact result, but it does not even go to zero for t —)■ oo. 
In that limit, one should have expected a resonating- 
tunneling bond at the site of the junction to trap an 
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FIG. 3. Comparison of the differential conductance for the 
simple junction calculated both consistently (or directly) and 
conventionally. Notice the unusual convention for the hori¬ 
zontal axis in order to highlight the t I/t duality of the 
problem. The vertical axis is in units of the single-channel 
quantum of conductance, Gq = jh, and G = "^Gq is the 
quantum limit for this problem. 


electron (for each spin) and thus block the passage of the 
current. A different way to describe it is by appealing 
to a tight-binding picture. The Hamiltonian for the two 
sites linked by t needs to be diagonalized first when t is 
the largest scale in the problem. One finds bonding and 
anti-bonding states that, in the t —> oo limit, will be al¬ 
ways occupied and always empty, respectively. The rest 
of the leads are relatively weakly coupled to these two 
states and not able to change their fillings and thus not 
able to produce a current. Instead of agreeing with this 
picture, the whole curve for the conventional TV char¬ 
acteristics resembles the result for a diode-like asymmet¬ 
ric junction with a non-Hermitian Hamiltonian [20], in 
which the formation of a resonating-tunneling bond is 
precluded by the model. This behavior alone could have 
been a clear indication that there are problems with the 
conventional way of calculating (even if one did not have 
a direct solution to compare with). 

The possible ramifications of our findings are many. 
A large number of the calculations done in the past, for 
example any problems sharing similarities with the one 
considered here (i.e., involving a junction, an impurity, or 
simply a boundary), will need to be reexamined critically. 
More generally, problems involving backward scattering 
or other types of nondiagonal interactions or processes 
need to be reconsidered for possible changes. Not all past 
results will be significantly affected though. For instance, 
on the one hand, the (weak-coupling) renormalization- 
group analysis of the effects of a “classical impurity” in 
a Tomonaga-Luttinger liquid [36] requires the knowledge 
of the impurity-potential beta function to leading order 
only, at which consistent and conventional calculations 
could be expected to (at least roughly) coincide; as our 


present results have shown it is the case if there are no 
interactions. There will be differences, but those would 
be expected in the finer details, probably appear at the 
next-leading order, and a calculation would be needed to 
determine them. On the other hand, the implications for 
the case of “quantum impurities” will be more dramatic. 
To put things in perspective, the good news is that we 
were able to provide a clear procedure, in the form of 
the n factors, to bosonize and debosonize a large class of 
models consistently. 

This paper was focused on the motivation and presen¬ 
tation of the formalistic details. In the future we will 
look at more involved examples of greater physical sig¬ 
nificance. We already started to reexamine some salient 
cases, and in the next paper we shall focus on the im¬ 
portant case of transport through quantum impurities in 
Fermi liquids [16]. 
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Appendix: Junction Thermodynamics 

The subtleties with the BdB-based mapping are quite 
generic and not restricted to nonequilibrium situations. 
Let us briefly compute the junction contribution to the 
free energy (or grand potential) and thus, indirectly, all 
thermodynamic quantities. We define the junction con¬ 
tribution in the same way as is done for impurity mod¬ 
els: as the difference between the full thermodynamic po¬ 
tentials with the junction closed and open, respectively. 
There is no voltage applied to the junction. 


1. Direct Calculation 

We start with the original (old) fermions and neglect 
the spin which will just give a factor of 2 at the end. Here 
we follow the procedure and notations as in Ref. 32. We 
will use, for convenience, a Nambu structure [otherwise 
we need to introduce sgn(w„) in the diagonal entries], 
but we do not need to use Keldysh and, instead, we will 
use Matsubara formalism and the following spinor basis: 

^ (Wn) = (i’h (W„) {-iOn) V'R {^n) V’r (“^n)) . 

Let us use again the definition 7 t = 2v-pt and use the 
standard result for the local inverse Green’s function for 
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the leads to write the local inverse Green’s function for With the same definitions, the local inverse Green’s func- 
the whole junction: tion for the junction is 


G-i (w„) 


—2ivp 


/I 0 zt 0 \ 

I 0 1 0 -it] 

lit* 0 1 0 

V 0 -it* 0 1 / 


(A.l) 


We compute the junction contribution to the thermody¬ 
namic potential via the standard method of “integrating 
over the coupling constant”: 


AO = - Oq = 





(A.2) 


Introducing the action determinant 


D (wn,^) = detGj ^ (w„) 

= \tf^^+2\t\^e+i 



we can use the formula 


AO = - 



E 

n^O 


d^D{0Jn,O 

D{uJn,0 


Since D{ojn,0 = -G(C) does not depend on frequency 
for the problem at hand, we factor out the divergent sum 
and indicate it as Sr=o = | J2n>o have 

= —Sr=0 In (^1 + 1^1 ^ 

^ -26r=o\n(l + \t\^) , (A.4) 

X spin \ / 


and we want to compare it with the result after the BdB- 
based transformations. 


2. Conventional Indirect Calculation 

Let us work in terms of the new fermions and adopt 
the following spinor basis: 


4 ' (w„) = (^tpi (w„) 2pj (-w„) -tjjsi (w„) V’L 


/I 0 -it* it*\ 

G-' = -2.W -f ’j . (A.5) 

\ it it* 0 1 / 

This time the action determinant reads 

0(w„,e) = 4|t|"e" + l, (A.6) 

and applying the same formulas we find 

AO = In (l + 4 |t|") . (A.7) 

But notice that if we “correct” the coupling constant we 
get 

AO -E^lnfl + ltl") . (A.8) 

We see that (i) the same “correction” as in the transport 
calculation is needed; (ii) we again lack an overall factor 
of 4, but (iii) the extra factor of (1 + |t| ) is not an is¬ 
sue this time (but notice that the logarithm would turn 
powers into factors). Because of the last point, a pertur¬ 
bative analysis is not effective to pinpoint the source of 
the discrepancies the way it is for transport calculations. 

Nota Bene'. Working in the consistent approach 
and using the (non number-eigenstates) half-filled basis 
(|0) ± |1)) /v^ in both the charge and spin sectors, if we 
have (he) = (hj) = l/-\/2 then the “correction” of the 
coupling constant reappears explicitly (but due to the 
use of CBCs). Moreover, tracing over the c and s sectors 
gives the missing factor of 4. One is thus able to recover 
the direct result with a calculation which does not differ 
much from the conventional one at the level of the local 
inverse Green’s function, but in an ad hoc way. 


3. Consistent Indirect Calculation 

Let us repeat the calculation but introducing yt —)■ 
Tto- = 7 thchs (recall we divided by 2 since we need to 
use CBCs). Using the appropriately modified result for 
the local inverse Green’s function for the junction, one 
finds the following action determinant: 

D (uj„, ^) = 1 -I- (2t*t^ + 2tlti) -I- (tj(t|)^ -I- t^(f*)^) 

= {\tfe+iy, ( a . 9 ) 

where the last expression is valid for either eigen- 
expectation-value of Ug. One thus recovers the same 
expression as in the direct calculation in the original- 
fermions language [cf. Eq. (A.3)]. All the ensuing results 
are thus identical. Notice the factor of 2 for spin will be 
contributed by tracing over eigenstates of h^, while on the 
charge sector only the (he) = 1 subspace contributes. 
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